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ABSTRACT 

The hard X-ray emission of active galactic nuclei (AGN) is believed to originate from the hot 
coronae above the cold accretion discs. The hard X-ray spectral index is found to be corre- 
lated with the Eddington ratio Lboi/^Edd, and the hard X-ray bolometric correction factor 
£boi/^x,2-iokcV increases with the Eddington ratio. The Compton reflection is also found to 
be correlated with the hard X-ray spectral index for Seyfert galaxies and X-ray binaries. These 
observational features provide very useful constraints on the accretion disc-corona model for 
AGN. We construct an accretion disc-corona model and calculate the spectra with different 
magnetic stress tensors in the cold discs, in which the corona is assumed to be heated by the 
reconnection of the magnetic fields generated by buoyancy instability in the cold accretion 
disc. Our calculations show that the magnetic stress tensor r rip — ap gas fails to explain all 
these observational features, while the disc-corona model with r rip = aptat always leads to 
constant Lboi/£x,2-iokcV independent of the Eddington ratio. The resulted spectra of the 
disc-corona systems with r rip = a^/p gas Ptot show that both the hard X-ray spectral index and 
the hard X-ray bolometric correction factor £boi/£x,2-iokcV increase with the Eddington 
ratio, which are qualitatively consistent with the observations. We find that the disc-corona 
model is unable to reproduce the observed very hard X-ray continuum emission from the 
sources accreting at low rates (e.g., Y ~ 1 for I/boi/^Edd ~ 0.01), which may imply the dif- 
ferent accretion mode in these low luminosity sources. We suggest that the disc-corona system 
transits to an advection-dominated accretion fiow+disc corona system at low accretion rates, 
which may be able to explain all the above-mentioned correlations. 

Key words: (galaxies:) quasars: general — accretion, accretion discs — black hole physics 



1 INTRODUCTION 

It is believed that active galactic nuclei (AGN) are powered by 
accretion of matter on to massive black holes, and the observed 
UV/optical emission of AGN is thought to be a thermal emission 
from the standard g eometrically thin, optically thick accretion disc s 



, optically I 

fe.g.. lShieldsll978| - |Malkan & Sargent! 198^ : 1 Sun & Malkari 19891) . 
On the other hand, the power-law hard X-ray spectra of AGN are 
most likely due to the inverse-Co mpton scattering of soft photons 
on a population of hot ele c trons dGaleev. Rosner & Vaiana 111 9791 ; 
Haard t & Maraschil Il99ll , 1 19931) . In the accretion disc-corona 
model, such soft photons are from the cold disc, a fraction 
of which are Compton scattered by the hot electrons in the 
corona above the cold disc to the hard X-ray energy band. The 
disc-corona m odel was extensively explore d in many previous 
works (e .g.. [Svensson & Zdziarskil 1 19941; lHaardt & Maraschil 



1991L 1 19931; iKawaguchi. Shimura & Mineshige! l200ll; 



Liu. Mineshige & Shibatal 20021 : lliu Mineshige & O hsuga 



magnetic fields generated in the cold disc are strongly buoyant, and 
a substantial fraction of magnetic energy is transported vertically 
to heat the corona above the disc with the reconnection of the 
fields (e.g.. |Pi Matted Il998l: iDi Matteo. Celotti & Fabiart Il999l : 
iMerloni & Fabian 200ll. |2002|) . It was found that the temperature 
of the hot electrons in the corona is roughly around 10 9 K, which 
can successfully reproduce a power-law hard X-ray spectrum as 
observed (e.g., Liu, Mineshige & Ohsuga 20031). 

The physical processes of turbulence triggered by the ampli- 
fied magnetic fields in the disc are very complicated, and quite 
unclear, though they are revealed to some exte nt with numerical 
magneto-hy drodynamical (MHD) simulations (Balb us~& Hawlevl 
Il99ll 11998). The so-called "a-prescri ption" is widely adopted i n 
most of the works on accretion discs jShakura & Sunvaevll 1973b . 
in which the magnetic stress tensor is assumed to be propor- 
tional to the total p ressure (p to t = Peas + Pr a d), gas pressure 
or -^/PgasPtot J Sakimoto & Coronitil Il98ll ; IStella & RosneJ 



20031) . In this disc-corona scenario, most gravitational energy is 



generated in the cold disc, probably t hrough turbulence pro duced 
by the magneto-rotational instability (Balb us & Hawlevll99ll) . The 



Pa 

1 19841 ; iTaam & Linl 1 1984 ). The magnetic stre s s tens or r rlp = 
aptot was suggested in IShakura & Sunvaevl dl973h . however, 
the disc is thermal unstabl e if the radiation pressu re domi- 
nates over the gas pressure dShakura & Sunvaevl 1 1976h (but also 
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see iHirose. Krolik & Blaesl l2008h . It was argued that the mag- 
netic field amplification is likely to be limited to that the mag- 
netic pressure is less than the gas pressure even in the in- 
ner radiation-pressure-dom inated regions, i.e., r r ^ = ap gas 
( Saki moto & Coronitill98lh. The disc w ith this stress tensor is ther- 
mal stable dSakimoto & Coroniti|[T98lh . The stress tensor r r ^ = 
oVPgasftot was initially suggested by iTaam &LirJdl984l) based 
on the idea that the viscosity is proportional to the gas pres- 
sure, but the size of turbulence should be limited by the disc 
scaleheight that is given in terms of the total pressure. This is 
also supported by the analysis on the local dynamical instabil- 
ities in magnetized, rad iation-pressure-supported accretion discs 
jBlaes &"So crates 2001), which also leads to p m = B 2 /8ty ~ 
Poy/PgasPtot, where p m is the magnetic pressure, and Bp is a con- 
stant of the order of unity (see iMerloni & F abian 2002j for a de- 
tailed discussion). iNavakshin. Rappaport & Melia I d2000l) argued 
that the accretion disc model with T rip = ap gas is too stable to ex- 
plain the unstable behaviour observed in GRS 1915+105, and they 
further suggested that it is possible to couple the radiation to the 
particles through collisions and thereby allow the radiation pressure 
to contribute to the stress tensor to some extent. As the complexity 
of the physics in the radiation-pressure-dominated fluids, any one 
of these stress tensors should only be regarded as a possible option 
in accretion discs. 

The hard X-ray observations on AGN may pro- 
vide useful clues on th e accr etion disc-corona models. 
IWang. Watarai & Mineshigel d2004l) compiled a sample of 
radio-quite AGN, and found a strong correlation between 
I/2-iokev/J'bQi and I/boi /I/Edd, which was conformed by 
Vasudevan & Fabian (2007) with a different AGN sample. This 
correlation was then used to constrain the disc-corona models with 
different magnetic stress tensors, and they found that the model 
with magnetic stress tenor r rip — ap gas is favored by the correla- 
tion of I/2-iokcv/iboi — iboi/iEdd- It was also found that the 
spectral index of hard X-ray continuum emission is correlated with 
the Eddington ratio dLu & Yul 1 1 9991 : IWang. Wata rai & Mineshige 
|2004 IShemmer et all 120061 l2008h . Izdziarski. Lubinski & Sniitnl 
(1999) found a strong correlation between the Compton reflec- 
tion and the hard X-ray spectral index for Seyfert galaxies and 
X-ray binaries. In the accretion disc-corona scenario, the hard 
X-ray emission originates from the Compton scattering of the 
soft photons by the hot electrons in the corona. Therefore, the 
correlations of the hard X-ray spectral index with the Edding- 
ton ratio/Compton reflection, together with the correlation of 
i2-iokcv/iboi — ^boi/iEdd, provide important constraints on 
the accretion disc-corona model. 

The parallel plane homogeneous corona is unable to produce 
an X-ray spectrum with T < 2 in 2-10 keV. The X-ray pho- 
tons radiated from the corona are reprocessed in the cold disc, 
and the reprocessed photons irradiate the corona, which cool the 
corona and lead to rather soft X-r ay spectra with T > 2. In the 
patchy corona model proposed by lHaardt. Maraschi & Ghisellinil 
( 1994), the corona appears as individual blobs above the cold disc. 
Most of the reprocessed photons do not enter the blob, so the 
cooling is significantly reduced and the blobs are hotter than the 
parallel plane homogeneous corona, which leads to harder X-ray 
spectra. This model can explain the observed hard X- ray spectra 
with r < 2 in some AGN (e.g.. Izdziarski et al.l 1996), however, 
it is still unable to explain the correlation between the Compton 
reflection and the har d X-ray spectral index for Seyfe rt galaxies 
and X-ray binaries (Zdzia rski. Lubinski & Smithll 19991) . which is 
also the case for the parallel plane homogeneous corona model. 



Izdziarski. Lubinski & Smith] dl999h suggested that a central hot 
plasma surrounded by a cold disc may explain the correlation of 
the Compton reflection with the hard X-ray spectral index. In this 
scenario, the cold disk is truncated at a certain radius d, within 
which the hot plasma may possibly c orrespond to an adv ection 
dominated accretion flow (ADAF) (e.g jNaravan & Yi|[i995r) . The 
radiation of the hot plasma is dominated by the inverse Compton 
scattering of the soft photons from the outer cold disc, while the 
cold disc is irradiated by the inner hot plasma. Thus, the Comp- 
ton reflection component decreases with increasing the inner ra- 
dius of the cold disc, which leads to less soft seed photons from 
the cold disc entering the hot plasma and then the harder X-ray 
spectrum from the hot plasma. The correlation between the Comp- 
ton reflection and the har d X-ray spectral index can be re pro- 
duced by this model (see Zdziarski, Lubinski & Smith 1999, for 
the de tails). An alternative model was proposed by Beloborodov 
( 1999) to explain the correlation between the Compton reflection 
and the X-ray spectral index. In this model, the hot plasma above 
the cold disc is assumed to move away from the cold disc at a mild 
relativistic velocity. Such an outflow reduces the downward flux, 
and then reduces both the reflection and reprocessing in the cold 
disc. The reduction of the reprocessing leads to less incident soft 
seed photons and cooling in the hot plasma above the cold disc , 
and in turn leads to harder X-ray spe ctrum l lBeloborodovlll999l : 
iMalzac, Beloborodov & Poutanenll200lh . 

In this work, we take these different magnetic stress tensors 
as candidates in our disc-corona model calculations, which could 
be tested with the X-ray observations of AGN. We summarize the 
disc-corona model in Sect. 2, and the numerical results of the model 
calculations are given in Sect. 3. In Sect. 4, we discuss the physical 
implications of the results. 



2 THE DISC-CORONA MODEL 

The gravitational power dissipated in unit surface area of the accre- 
tion disc is given by 



on 



1 



1/2' 



(1) 



where M is the mass accretion rate of the dis c, Qk(R) is the Ke- 
pleria n velocity at radius R, and R ln = 3Rs dShakura & Sunvaevl 
1 19731) . The Schwarzschild radius R s = 2GM bh /c 2 , where M bh is 
the black hole mass. The corona is assumed to be heated by the re- 
connection of the magnetic fields generated by buoyancy instability 

in the disc. 

The power dissipated in the corona is dDiMatteoll 19981) 



r»- - - B 

/cor — Pm'fp — n ^p, 
07T 



(2) 



where p m is the magnetic pressure in the disc, and v p is the veloc- 
ity of the magnetic flux transported vertically in the disc. The ris- 
ing speed v p is assumed to be proportional to their internal Alfven 
velocity, i.e., v p = bv\, in which b is of the order of unity for 
extremely evacuated magnetic tubes. 

The soft photons from the disc are Compton scattered by the 
hot electrons in the corona to X-ray bands, and about half of the 
scattered photons are intercepted by the disc. The reflection albedo 
a is relatively low, a ~ 0.1 — 0.2, and most of the incident pho- 
tons from the corona are re-radia ted as blackbody radiation (e.g., 
Izdziarski. Lubinski & Smithll 19991) . Thus, the energy equation for 
the cold disc is 
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+ -(l-a)Q+ r = 



4aT^ 
3r 



(3) 



where Tdi sc is the effective temperature in the mid-plane of the disc, 
and t = r cs + Tg is the optical depth in vertical direction of the 
disc. In this work, we adopt a = 0.15 in all our calculations. 

As the detailed physics for generating magnetic fields in the 
accretion disc is still quite unclear, we adopt different magnetic 
stress tensors as: 



(4) 



OyPgasftot, 



in our model calculations, respectively. We summarize the equa- 
tions describing the disc as follows: 

The continuity equation of the disc is 



4irRH d (R)p(R) Vl{ (R) = M, 



(5) 



where Ha(R) is the half thickness of the disc, p(R) is the mean 
density of the disc, and v-r(R) is the radial velocity of the accretion 
flow at radius R. 

The equation of state for the gas in the disc is 



Ptot = 



+ Prad = 



+ g a Tdi sc , 



(6) 



where p — (1/tM + l/fi e )~ , pi = 1.23 and p c — 1.14 are 
adopted corresponding to the plasma consisting of 3/4 hydrogen 
and 1/4 helium. The vertical hydrodynamical equilibrium requires 
c b /Qk, where the sound speed c s 



H, 



a/ (Ptot +Pm)/P- 



The angular momentum equation for the disc is 
i?i^ 1/2 l 



MQk(R) 



1 - 



R 



= A-kH a t tv , 



(7) 



where the magnetic stress tensor is given by equation l[4}. 

Solving equations l[T}-(|7) numerically, the structure of the disc 
and the power dissipated in the corona Qtor can be derived as func- 
tions of radius R. The ratio of the power dissipated in the corona to 
the total for such a disc-corona system is available: 



</> = 



/ Q+ OI 2-nRAR 



Pcor 



J Q+ iss; 27rRdR' 
The equation of state for the hot gas in the corona is 



fnm p 



p c m p 



+ Pcor 



(8) 



(9) 



where T and T c are the temperatures of the ions and elec- 
trons in the two-temperature corona, and the magnetic pressure 
Pcor.m = -Bcor/87r. In this work, the magnetic fields are assumed 
to be equipartition with the gas pressure in the corona. The ver- 
tical hydrodynamical equilibrium in the corona requires H COI = 
Ccor, s /^K, where the sound speed c co 



- ^/pcor / pcor • 

The energy equation describing the two-temperature corona is 



Qc 



FZ 



(10) 



fsyn + F htem + F Com P is the cooling rate in unit sur- 



where F cor 

face area of the corona and Qcor is the energy transfer rate from the 
ions to the electrons in the corona via Coulomb collisions, which 
is given bv lStepnev & Guilberj fl983h . The fraction of the energy 
directly heat the electrons 5 can be as high as ~ 0.5 by mag- 
netic reconnection, if the magnetic fields i n the plasma are strong 
jBisnovatvi- Kogan & Lovelacefr9 97. 2000). In this work, we adopt 
8 = 0.5 in our calculations. For the plasma consisting of different 



elements, the C oulomb interact ion between the electrons and ions 
is given by (see[zdziarski 1 998h 



QcOT 



1.5 



E 



m e 2 



kTi - kT e 



2(e c + Q Z ) 2 + i„ /e c + e 



e c + &i 



K 2 {i/e e )K 2 (i/ez) 
e c + ez 



In A 



,(ii) 



where Az is the mass number of the Zth element, In A = 20, 
6 C = kT c /m c c 2 , Q z = ®i/Az, and 9 ; = kTi/m p c 2 are the di- 
mensionless temperatures. Although the cooling rate of the corona 
is dominated by the Compton scattering of the soft incident pho- 
tons from the disc, we include the synchrotron, bremsstrahlung and 
Compton emission in our calculations. The cooling terms F S y n and 
F brem are the functions of electron number density, temperature, 
and the ma gnetic field strength o f the gas in the corona, which are 
taken from lNaravan & Yil j 19951) . 

Assuming the corona to be a parallel plane, the mean proba- 
bility of the soft photons injected from the cold disc experiencing 
the first-order scattering in the corona is 



Pi 



(1- 



-To / COS(?\ 



cos 6d cos 8, 



(12) 



where the constant specific intensity of the soft photons from the 
disc is assumed, 8 is the angle of the motion of the soft photons with 
respect to the vertical direction of the disc and To = aTn e H COT is 
the optical depth of the corona for electron scattering in the vertical 
direction. For simplicity, we assume the first-order scattered pho- 
tons are radiated uniformly throughout the vertical direction of the 
corona. The mean probability for these first-order scattered photons 
experiencing the second-order scattering can be estimated as 



-£tq / cos G 



+1- 



-(l-€)Wcos0 



]d cos 0,(13) 



where £ = z/H COI , and the first-order scattered photons are as- 
sumed to be radiated isotropically. The probability of the scattered 
photons experiencing next higher order scattering approximates to 
P2, so we simply ad opt P n = P2 for n > 2 . Thus, using the 
method suggested by Coppi & Blandford (1990), we can calculate 
the Comptonized spectrum with equations | |12I > and J 1 3b . if the den- 
sity, temperature of electrons, and the incident spectrum of the disc 
are known. Integrating the derived Comptonized spectrum over the 
frequency, the cooling of the electrons in the corona due to the 
Compton scattering Fq is available. 

The disc structure and the power dissipated in the corona 
can be calculated as functions of R, provided the black hole 
mass Mbh, mass accretion rate M and the viscosity parameter a 
are specified. The scaleheight of the corona is a function of Ti, 
T e , and n e , on the assumption of static hydrodynamical equilib- 
rium in the vertical direction. There are two equations describing 
the energy equilibrium in the corona at R (equation |10t , which 
contain three physical quantities of the corona: Ti, T c , and n a . 
Liu, Mineshige & Ohsuga ( 2003) calculated the vertical structure 
of the two-temperature corona above a thin disc based on the evapo- 
ration mechanism, and their results show that the temperature of the 
ions in the corona is always in the range of ~ 0.2 — 0.3 virial tem- 
perature (the virial temperature is defined as T v i r = GMm p /kR 
in their work). In this work, we adopt a slightl y different defini- 
tion: T v ir = GMm p /3kR, as that adopted by Merl oni & Fabianl 
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1 I2002I) . To avoid the complexity of calculating the vertical struc- 
ture o f the corona, we simply adopt T j = 0.9T v i r [equivalent to 
0.3 in lLiu. Mineshige & Ohsugal d2003l) 's work] in our model cal- 
culations. Thus, the structure of the corona, i.e., the temperature 
and density of the electrons, can be derived, and the spectra of the 
disc-corona system are available based on the derived disc-corona 
structure with different magnetic stress tensors given in equation 



3 RESULTS 

We calculate the disc-corona structure as described in section 2. 
The black hole mass Mbh = 1O 8 M0 is adopted in all our cal- 
culations, because the main features of the hard X-ray spectra are 
almost independent of Mbh for massive black holes in AGN. In 
Fig.Q] we plot the ratios of the power radiated in the corona to the 
total power L C ar/£bol as functions of accretion rate rh predicted 
by the models with different magnetic stress tensors. The dimen- 
sionless accretion rate rh is defined as rh = M /Msaa, where 
Msdd = Lvdd/Vesc 2 , and a conventional radiative efficiency 
VaB = 0.1 is adopted. The model with r rip = ap to t always leads to 
constant ratios L cor /Lboi independent of accretion rate rh, while 
the disc-corona model calculations with T rip = ap gas show that 
Lcoi/ ibol becomes extremely small for high accretion rates (e.g., 
L COI /L ho i ~ 0.01 form ~ 1). The model with r r ^ = Q^/PgasPtot 
shows that the ratios L cor /Lb \ ~ 0.3 — 0.6 at rh = 0.01, while 
^cor/iboi ~ 0.05 — 0.1 at rh = 1, for different values of a. 

The temperature and optical depth for Compton scattering of 
the hot electrons in the vertical direction of the corona are plotted 
in Fig. [2] The electron temperatures are in the range of ~ 5 x 
10 8 — 3 x 10 9 K for different values of rh. The temperature of 
the hot electrons in the corona tends to decrease with accretion rate 
rh. When the accretion rate rh is as high as ~ 0.5, the electron 
temperature of the corona decreases to ~ 5 x 10 s K. 

In Fig. [3] we plot the spectra of the disc-corona systems calcu- 
lated with different magnetic stress tensors. The hard X-ray emis- 
sion indeed exhibits a power-law feature for all the models. The 
photon spectral indices and the ratio of the bolometric luminos- 
ity to the X-ray luminosity in 2-10 keV as functions of accretion 
rate for different magnetic stress models are plotted in Fig. [4] The 
photon indices T do not change much with accretion rate rh for 



the disc-corona models with either r r 



aptot or T r , 



aPgas, 



while the hard X-ray spectral index T increases significantly with 
accretion rate rh for the model with r rip = a^/p gaB p tot . 

We compare the spectra of the accretion disc/corona systems 
with different ion temperature 7] adopted in Fig. [5] It is found 
that the X-ray spectra change little with the value of Ti, provided 
all other parameters are fixed. In Fig. [6] we plot the Compton y- 
parameter (y — 4kT c To/m e c 2 ) varying with 7], and find that it is 
insensitive to T\ . 




Figure 1. The ratios L cor /Lbol as functions of accretion rate rh predicted 
by the models with different magnetic stress tensors. The red lines represent 
the results calculated with the magnetic stress tensor T rip = aptot (green 
lines: r rlp = crp gas ; blue lines: r rip = ct^PgasPtot)- The different line 
types represent the different values of viscosity parameter a adopted (solid 
lines: a = 0.2; dashed lines: a = 0.3; dash-dotted lines: a = 0.5 and 
dotted lines: a = 1). 
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4 DISCUSSION 

It is believed that the power generated in the disc is transported 
vertically with the buoyancy of the magnetic fields in the disc, 
and the fraction of the power dissipated in the corona to the to- 
tal is mainly regulated by the magnetic fields (e.g., Di Matted 
1998l:lMerloni & Fabianl2002l ; IWang, Watarai & MineshigekOoUT 
Vasudevan & Fabianl (2007) found that the fraction L COI /Lhoi is 
about 0.5 for the sources with low Eddington ratio Lbol/^Edd ~ 



Figure 2. The upper panel: the temperature of the electrons in the corona as 
functions of disc radius R. The colours represent the models with different 
magnetic stress tensors, which are the same as those in Fig. [T] The solid 
lines represent the results calculated for rh = 0.5, while the dashed and 
dash-dotted lines are for rh = 0.05 and 0.005, respectively. The viscosity 
parameter a = 0.5 is adopted for all the model calculations. In the lower 
panel, we plot the optical depth for the Compton scattering of the electrons 
in the vertical direction of the corona. 
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10 10 
v(Hz) 

Figure 3. The spectra of the disc-corona systems with different magnetic 
stress tensors (a: r rip = aptot; b: r rlp = ap gas ; c: r rLp = a^/PgasPtot). 
In all the calculations, a = 0.5 is adopted. The dotted lines represent the 
spectra of the cold discs, while the dashed and dash-dotted lines are for 
the Compton and synchrotron+bremsstrahlung radiation respectively. The 
accretion rates: rh = 0.5(red), 0.05(green) and 0.005(blue) are adopted 
respectively(from up to down). The two black dotted lines represent 2 and 
10 keV, respectively. 



0.01, while it decreases to 0.1 for Lboi/^Edd ~ 1, for a sample of 
AGN. From Fig. Q] we find that the model with r rip = ap to t al- 
ways predict a constant L COI /Lboi, which seems to be inconsistent 
with the observation. The models with other two magnetic stress 
tensors can roughly reproduce the trend of Z/ cor /Z/Edd decreasing 
with accretion rate m, though the model with r rip — ap ga s always 
underpredicts the L cor /I/boi at high Eddington ratio end, even if 
q = 1 is adopted. It seems that the model with r rv = a^/p g asPtot 
is better than the other two models. 

The ratio L C or/£boi decreases with accretion rate rh, which 
means the fraction of the power radiated from the cold disc in- 
creases with rh. The cooling mechanism of the hot corona is dom- 
inated by the inverse Compton scattering of the soft photons from 
the cold disc by the hot electrons in the corona. There are more soft 
photons supplied by the cold disc for high-m cases, and the frac- 
tion of the power radiated from the corona decreases with rh. Thus, 
it can be easily understood that the temperature of the hot electrons 
in the corona should decrease with accretion rate rh (see Fig. [2j. 
It is found that the temperature of the electrons in the corona de- 
creases to ~ 



— 4 x 10 K for the models with r rv = ap gas 
or Trip = Q^/PgasPtot when rh — 0.5, which implies that the hot 
corona is nearly to be suppressed when rh is as high as > 0.5. 

We calculate the spectra of the accretion disc-corona systems, 
and find that their hard X-ray continuum emission indeed exhibits 
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Figure 4. The upper panel: the photon spectral indices as functions of ac- 
cretion rate rh for different magnetic stress models (red: T rlp = aptot; 
green: r rtp = ap ga s; blue: T rtp = ce^/PgasPtot). The squares represent 
the results calculated with a = 0.5, while the triangles are for a = 0.3. 
The X-ray bolometric correction factors Lbol/^x 2— lOkeV as functions of 
accretion rate rh for different models are plotted in the lower panel. 




Figure 5. The spectra of the disc corona systems calculated for rh = 0.01, 
and 0.1, respectively (Trip = ce^/pgasptot is adopted in the calculations). 
The red lines represent the spectra for Tj = 0.9T v ; r , while the green 
lines are for Tj = 0.3T v ; r . The dotted lines represent the spectra of the 
cold discs, while the dashed and dash-dotted lines are for the Compton and 
synchrotron+bremsstrahlung radiation respectively. 



a power-law spectrum (see Fig. f5). In Fig. 4, we plot the photon 
spectral indices F in 2-10 keV and the X-ray correction factors 
£boi/£x,2-iokcv as functions of accretion rate rh with different 
magnetic stress tensors. We find that the photon spectral index V 
predicted by the model with Trip = aptot is within ~ 2 — 2.2 for 
different values of rh, and it increases slightly with rh. The X-ray 
correction factor Lbol/£x,2-iokaV remains constant with accre- 
tion rate rh for this model. The photon spectral index T increases 
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Figure 6. The Compton {/-parameter versus the ion temperature Tj in the 
corona at radius R. with different accretion rates m = O.Ol(red), and 
O.l(green). The solid lines are calculated with R = 10/?g, while the dashed 
lines are for R = 40i?s . 



with m while rh < 0.02, and it then decreases with rh, for the 
model with r rip — ap gas . We find that the bremsstrahlung emis- 
sion contributes some to the hard X-ray energy band when rh is 
high, and the photon spectral index T is slightly affected by the 
bremsstrahlung emission (see the middle panel of Fig. 3). We find 
that the X-ray correction factor I/boi/ix,2-iokeV f° r this model 
becomes extremely high for high-m cases, which seems to be in- 
consistent with that derive d from the observations (see Fig. 12 in 
IVasudevan & Fabian| [2007). The X-ray spectra of the disc-corona 
systems with r r ^ = a^/pgasptot show that both the photon spectral 
index V and the X-ray correction factor -Lboi/ix,2-iokcV increase 
with accretion rate m, which are qualitatively consistent with the 
observations, except for the low rh end. The temperature of the 
electrons in the corona decreases with rh (see Fig. |2j, and the in- 
verse Compton scattered X-ray spectrum therefore becomes softer 
for a higher rh. Similarly, anticorrelations between the X-ray spec- 
tral i ndices and X-ray fluxes were found i n some X-ray binarie s 
(e.g., lYu et alj|2003l; IZdziarski et ai]|2004 iDel Santo et ai]|2008l) . 
which seem to be roughly consistent with our present disc corona 
model. 

The photon spectral index V is found to increase with the Ed- 
dington ratio, and it could b e as low as ~ 1 when rh ~ 0.01 (see 
Shemmer et al. 2006, 2008h . The hard X-ray photon spectral index 
r ~ 2 derived from our model calculations with any magnetic 
stress tensor when rh ~ 0.01, which are significan tly higher than 
the observed V ~ 1 (e.g. lShemmer et al.l200ril2008l). Although the 
atchy corona model proposed by lHaardt, Maraschi & Ghisellinil 
19941) can produce the observed hard X-ray spectra with 1? ~ 1, 



this patchy corona model is unable to explain the correlation be- 
tween the Compton reflection and the hard X-ray spectral index. In 
an alternative model, the hot plasma above the cold disc is assumed 
to move away from the cold disc at a mild relativistic velocity 
dBeloborodovlll999h . which reduces both the reflection and repro- 
cessing in the cold disc. This naturally leads to harder X-ray spec- 
trum (see the discussion in Sect. 1). However, the detailed physics 
for producing such mildly relativistic outflows is still unclear. 

It was argued that the central engines in these low-luminosity 
sources with very hard X-ray emission may be different from their 



high-luminosity counterparts, i.e., the advection dominated accre- 
tion flows (ADA Fs) may be present in these low-luminosity sources 
dLu & Yu1ll999h . The soft incident photons for Comptonization are 
mainly due to the synchrotron+bremsstrahlung emission in the 
ADAFs. The energy density of the soft photons in the ADAF is 
much lower than that in the disc-corona model, which leads to in- 
efficient cool ing and relatively hig her electron temperature in the 
ADAFs (e.g jNaravan & Yilfl995l) . Thus, the X-ray spectra of the 
ADAFs can be much harder than those of the disc-corona sys- 
tems. The very hard X-ray spectra observed in these low lumi- 
nosity sources can be well modelled with the ADAFs accr eting at 
low rates (see e.g., lOuataertet all 19991: IXu & Cad 2008, for the 
spectral modeling for the hard X-ray emission from the low lu- 
minosity AGN with ADAFs). Our present model calculations are 
limited to the disc-corona model, in which the cold discs extend 
to the marginal stable orbits. The model of a disc-corona connect- 
ing with inner ADAF may resolve the photon index discrepancy at 
the low rh end. The geometry of this ADAF+disc/corona scenario 
is quite similar to the hot plasma+ cold disc model proposed by 
IZdziarski. Lubinski & Smith] d 19991) . When the accretion rate de- 
creases to a critical value m cr it, the inner cold disc may be trun- 
cated at radius R tI and it transits to an ADAF within this radius. 
The transition radius R tI may be close to the mar ginal stable or- 
bits soon after the accretion mode transition (e .g., lOuataert et al.l 
1999: lYuan & Naravanl 120041 ; IXu & Caol koOS), and then it may 
increase with decreasing accretion rate rh (e.g ., Liu et al.l 19991; 



i Rozanska & Czemvl200(j ; ISpruit & Deufell2002l ; [Yuan & Naravanl 
2004). The X-ray spectrum of such an ADAF+disc/corona system 
consists of emission from the inner ADAF and outer corona. The 
ratio of the X-ray emission from the ADAF to that from the corona 
increases with decreasing m, and therefore the photon spectral in- 
dex r may decrease smoothly with decreasing rh provided the ini- 
tial truncated radius of the cold disc is not very large compared with 
the ma rginal stable orbits. Similar to IZdziarski. Lubinski & Smith! 
( 1999)'s model, the correlation between the Compton reflection 
and the hard X-ray spectral index can also be naturally explained 
by this ADAF+disc/corona model, if the truncated radius R tr in- 
creases with decreasing accretion rate rh. This model can also 
succ essfully explain the spectral be haviours of X-ray binaries 
(see lDone, Gierlinski & Kubotall2007l . for a review and references 
therein). The detailed calculations on such ADAF+corona systems 
will be reported in our future work. 

In this work, we simply assume 71 = 0.9T v i r , motivated 
by the previous work on the disc-corona model calculation s 
dLiu. Mineshige & Shibatal200ilLiuT Mineshige & Ohsugal2003l) . 
We also check how 71 may affect our results by tuning the value 
of 71, and find that the X-ray spectra of the disc/corona systems 
change very little if all other disc parameters are fixed (see Fig.[5jl. 
The cooling of the corona is dominated by the inverse Compton 
radiation, which is roughly proportional to Compton y-parameter. 
Thus, it is not surprising that our calculations show the Compton y- 
parameter varying little with 7] if all other parameters are fixed (see 
equation [TO] and Fig.[6}. The thickness of the corona is mainly reg- 
ulated by the ion temperature 71, and therefore the electron density 
decreases with increasing 71, which leads to significant difference 
in synchrotron/bremsstrahlung spectra for different 71 (see Fig. [5}. 
We find that the resulted rh — T relations are almost not changed if 
a different value of 71 is adopted. Unlike the ADAFs, almost all the 
power dissipated in the hot corona with magnetic reconnection is 
radiated away locally. This means the radiated power in the corona 
is independent of the value S, and the temperature and density of 
the electrons in the corona are almost insensitive with the value of 
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5. We also perform the same model calculations for different black 
hole masses (e.g., Mbh = 1O 9 M0). It is found that our results are 
almost independent of Mbh for massive black holes. 

In all our calculations, we assume the magnetic fields to be 
equipartitioned with the gas pressure in the corona. As the cool- 
ing of the corona is dominated by the inverse Comptonization 
of the soft photons from the cold disc, the structure of the disc- 
corona and its spectrum (except in the radio wavebands) is al- 
most i ndependent of t he ma gnetic field strength in the corona. Re- 
cently, Laor &BehaJ [2008) found there is a strong correlation be- 
tween the radio luminosity (Lr) and X-ray luminosity (Lx) with 
Lr, ~ 10~ 5 Z/x, for the radio quiet Palomar-Green (PG) quasar 
sample. The spectra of our disc-corona model show that the radio 
emission is correlated with the X-ray emission, which is roughly 
consistent with the c orrelation between Lr and Lx discovered by 
lLaor&Beha3 d2008h . 
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